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ABSTRACT 

Millisecond pulsars (MSPs) have been firmly established as a class of 7-ray emitters 
via the detection of pulsations above 0.1 GeV from eight MSPs by the Fermi Large Area 
Telescope (LAT). Using thirteen months of LAT data significant 7-ray pulsations at the 
radio period have been detected from the MSP PSR J0034— 0534, making it the ninth 
clear MSP detection by the LAT. The 7-ray light curve shows two peaks separated by 
0.274±0.015 in phase which are very nearly aligned with the radio peaks, a phenomenon 
seen only in the Crab pulsar until now. The >0.1 GeV spectrum of this pulsar is well 
fit by an exponentially cutoff power law with a cutoff energy of 1.8±0.6±0.1 GeV and a 
photon index of 1.5ib0.2±0.1, first errors are statistical and second are systematic. The 
near-alignment of the radio and 7-ray peaks strongly suggests that the radio and 7-ray 
emission regions are co-located and both are the result of caustic formation. 

Subject headings: Facilities: Fermi () Pulsars: gamma-ray, individual (PSR J0034— 0534) 



1. INTRODUCTION 



Forty new pulsars have recently been observed to pulse in high-energy (HE), >0.1 GeV, 7- 
rays by t he Large Area Telesc ope (LAT) aboard the Fermi Gamma-ray Space Telescope (formerly 
GLAST) (lAbdo et al.ll20010al ). Among these new detections are eight millisecond pul sars (MSPs 



(lAbdo et al.ll2009al ). MSPs are thought to be older, recycled pulsars in binary systems (jAlpar et al 



a theory which has been supported by the discovery of millisecond pulsations from accreting 
low-mass X-ray binaries (LMXBs) (e.g. Wijnands & van der Klis 1998) and more recently b y 
obser v ations of an MSP transitioning from the LMXB to the pulsar phase (lArchibald et al.ll2009l ). 
Usovl (jl983l ) first attempted to show that MSPs should be 7-ray emitters and could, possibly, 
have an even greater HE luminosity than the pulsar in the Cr ab nebula. Data from the Energetic 
Gamma-Ray Experiment Telescope (EGR ET) dThompsonlboOS ) were searched for both point source 
and pulsed emission from several MSPs (|Fierro et al.lll995l ) and, while n one were dete c ted, u pper 
limits were set which put useful constraints on MSP emission models. iKuiper et al.l (|2000l ) did 
report a mar ginal (^ 4cr) pulsed detection of the MSP PSR J0218+4232, which has been confirmed 
with Fermi (jAbdo et al.ll2009al ). HE 7-ray emission has also been detected from the vicinity of the 
globular cluster 47 Tucanae with a 7-ray spectrum coi isistent, at the 95% confidence level, with the 
superposition of between seven and sixty- two MSPs (lAbdo et al.ll2009bl ). The Italian observatory 



^^Consorzio Intemniversitario per la Fisica Spaziale (CIFS), 1-10133 Torino, Italy 
North- West University, Potchefstroom Campus, Potchefstroom 2520, South Africa 



59 



Dipartimento di Fisica, Universita di Roma "Tor Vergata", 1-00133 Roma, Italy 



■'School of Pure and Applied Natural Sciences, University of Kalmar, SE-391 82 Kalmar, Sweden 
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AGILE has also reported a 4.2cr d etection of 7-ray pulsations from PSR B1821— 24 in the globular 
cluster M28 dPellizzoni et al.1 bood ') . 



Fermi began nominal sky-survey observations on 2008 August 4, viewing the entire sky every 
two orbits (~3 hours). The main instrument on Fermi is the Large Area Telescope f LAT), which is 



a pai r-conversion telescope sensitive to 7-rays with energies from 0.02 to >300 GeV (jAtwood et al, 
20091 ). The LAT has a 2.4 sr field of view, a peak effective area of ~ 8000 cm^ above 1 GeV on axis 



and a 68% containment radius of 0.6° at 1 GeV for events converting in the front section of the 
LAT. The LAT timing is derived from a GPS clock on the spacecraft, and events are time-stamped 
to an accuracy better than 1 (jAbdo et al.ll2009d ). Using approximately thirteen months of LAT 
data, we have discovered 7-ray emission from PSR J0034— 0534, which thus becomes the ninth 
7-ray MSP detected with by Fermi LAT. 



PSR J0034-0534 



PSR J0034— 0534 was discovered by iBailes et al.l (j 19941 ) in a survey of the southern sky with 
the Parkes radio telescope. The pulsar is in a binary sytem with a 1.6 d orbital period. The initial 
observations measured a spin period (P) of 1.877 ms, which makes this the fastest 7-ray MSP 
yet detected with the LAT, and a period derivative (P) of 6.7 x 10^^^. Table 1 lists some of the 
measured and derived timing parameters for PSR J0034— 0534 using the radio data described in 
§3.1. The derived quantities in Table 1 are corre cted for the Shklov skii effect (jShklovskiil Il970l ) 



and assume a moment of inertia I = 10^^ g cm^. lAbdo et al.l ( 2009al ) observed that 7-ray MSPs 
had similar spectral properties and, due to their short periods, comparable m agnetic field s trength s 
at the light cylinder (Blc) to those of younger 7-ray pulsars. Additionally, lAbdo et al.l (l20010al ) 
noted a weak correlation of cutoff energy with Blc for all forty-six 7-ray pulsars detected by the 
LAT in the first six months. Of the 7-ray MSPs detected so far only PSR J0218+4232 has a higher 
value of Blc than PSR J0034— 0534 while most of the younger 7-ray pulsars detected with Fermi 
have lower values of Blc • 

Large values of Blc have been linked to giant pulses (GP) in the radio (Cognard et al. 1996 
and Knight et al. 2005) and the relatively large value for PSR J0034— 0534 prompted several GP 
searches (Romani & Johnston 2001, McLaughlin &: Cordes 2003, and Knight et al. 2005) though 
none were detected. Hubble Space Telescope observations revealed an optical co unterpart to the 



binar y companion of PSR J0034— 0534 consistent with a white dwarf hypothesis (jLundgren et al 
19961 ). Infrared observations of PSR J0034— 0534 only put upper limits on the flux density of 
any surronding debris disk (Greaves & Holland 2000 and Lazio &: Fischer 2004). XMM-Newton 
observations revealed a low significance (< 3a) X-ray source 0.2" from the radio position of PSR 
J0034— 05 34, but no pu lsed signal was detected from the source and a firm identification could not 
be made (|Zavlinl bood ). The EGRET 3a point source flux upper limit, 0.1 to 10 GeV, for PSR 
J0034-0534 of 15.2 x 10~^ cm"^ s~^ (jPierro et al.lll995l ) is an order of magnitude above the LAT 
flux in S4.2. 
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OBSERVATIONS 



3.1. Radio Timing 



The timing solution used for PSR J 0034— 0534 has been derived from observations carried out 
at the NanQay radio telescope in Prance (ITheureau et al.ll2005r) and t he Westerbork Synthesis Radio 
Telescope (WSRT) in the Netherlands ([Karuppusamv et al.l l2008l ). using 170 Times Of Arrival 
(TOAs) recorded between 2005 October 26 and 2009 August 24. Among the 170, 139 were recorded 
at the WSRT at frequencies b etween 314 MHz and 376 MHz with a mean uncertainty of 2.5 /is 
and a bandwidth of 10 MHz (jStappers et al.ll20ld ). The remaining TOAs were recorded at the 
Nangay radio telescope at 1398 MHz with a mean uncertainty of 3.2 u s with a bandwidth of 64 
MHz before 2008 June 13 and 128 MHz thereafter (jCognard et al.ll2009l ). This multiple-frequency 
dataset tightly constrains the dispersion measure (DM), which is crucial for profile comparisons 
at different wavele ngths. The ephemeris was derived using the TEMPO^ pulsar timing package 
(jHobbs et al.ll2006l ). fitting for the rotation frequency and its first derivative while accounting for 
the binary motion of the millisecond pulsar. The combined timing solution gives a post-fit rms of 
7.1 ^s and a DM measurement of 13.76517 it 0.00004 cm~'^ pc, with no indication of variation with 
time. The uncertainty in the dispersion measure leads to an uncertainty of less than 2 us in the 
extrapolation of 300 MHz arrival times to infinite frequency, negligible for the low-s tatistics gamma- 
ray li ghtcurve of PSR J0034— 0534. Using the NE2001 electron density model of ICordes fc: Lazio 
( 200^ 11 and the measured DM gives a distance of d = 0.53 =fc 0.21 kp c, assuming 40% uncertainty 
due to fluctuations in the free electron density (jBrisken et al.l 12002) • With this dataset it was 
not possible to measure the timing parallax and extract a more accurate distance measurement. 
However, for a pulsar at a distance of ~0.5 k pc and near the ecliptic plane the contribution of 
parallax to timing residuals is less than 2.4 fj,s (jLorimer &: Kramerl |2004| ) and thus well below the 
precision of this ephemeris. This timing solution will be made available through the Fermi Science 
Support CenteJil. 



3.2. LAT Data Selection 

Analysis of LAT data was done using the Fermi Science Tool^ (STs) v9rl5p2. The Fermi ST 
gtselect was used to select events which had reconstructed sky directions within 10° of the radio 
position of PSR J0034— 0534, energies from 0.1 to 100 GeV, zenith angles < 105°, and spanning 
2008 August 4 to 2009 September 10. Events were required to belong to the "Diffuse" class of 



*http: / /tempo2. sourceforge.net / 

^ http://rsd-www.nrl. navy.mil/7213/lazio/ne_model/ 

*littp://fermi. gsfc.nasa.gov/ssc/data/access/lat/ephems/ 

^ http://fermi.gsfc. nasa.gov/ssc/data/analysis/scitools/overview. html 
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events as defined under the P6_V3 instrunient re sponse functions (IRFs), those with the highest 
probability of being photons (jAtwood et al.ll2009l ). Additionally, the Fermi ST gtmktime was used 
to exclude times when the rocking angle of the instrument exceeded 52° and when the Earth's limb 
infringed upon the 10° region of interest. The events were then phase folded with the radio timing 
solution using the Fermi plug-ir0 now provided with the TEMP02 software. 



4. RESULTS 
4.1. Light Curve 

Events found within 0.8° of PSR J0034— 0534 were selected from the LAT data described in 



3.2 and tested for periodicity with the Fermi ST gtptest. The result was an H-test (jde Jager et al 



19891 ) value of 47.7 with a chance probability of 6.8 x 10~ , corresponding to a pulsed detection 
significance of 5.4cj. Figure 1 shows the folded light curve of these events over two rotation periods 
for events above 0.1 GeV and 1 GeV as well as the 300 MHz WSRT and 1.4 GHz Nancay radio 
profiles. There is a small contribution to the pulse width in the radio profiles due to scattering; 
however, this contribution is much less than the bin width used in Figure 1. The 7-ray light curve 
of PSR J0034— 0534 shows two peaks which are very nearly aligned with the radio peaks. This 
near-alignment of the radio and 7-ray peaks is very reminiscent of what is seen in the Crab pulsar 
(e.g. Abdo et al. 2010b) which is nearly aligned in radio, optical. X-ray, and 7-rays (GeV and TeV). 
The >0.1 GeV light curve in Figure 1 was fit with two Lorentzians plus a constant offset, fixed at 
the value of the background estimate shown, which gave a reduced value of ~ 1.4 indicating 
good agreement with the data. Table 2 lists the peak positions {(pi), full width half maximum 
(FWHMi), radio to 7-ray phase lags {5i), and peak separation (A) values. The phase lag values in 
Table 2 were calculated by assuming the first radio peak to be at phase and estimating the second 
at phase 0.258, using the 324 MHz radio profile. The phase lag in the second peak is statistically 
consistent with while the first peak has a significant, but small, offset from the radio. Fitting the 
7-ray light curve with asymmetric Lorentzians did not improve the fit, though with more data the 
second peak may show significant asymmetry. Using only events >1.4 GeV gives a pulsed detection 
of 3a while using events >2 GeV gives a pulsed detection of only 1.6a. This indicates that there is 
significant evidence for emission up to almost 2 GeV from PSR J0034— 0534. 



4.2. Spectrum 

An unbinned maximum likelihood method (Cash 1979 and Mattox et al. 1996), using the 
pyLikelihood python module included with the Fermi STs, was used to fit the region around PSR 
J0034— 0534. All point sources found above the background with test statistic >25 in a preliminary 



http://fermi.gsfc.nasa.g0v/ssc/data/analysis/scitools/pulsar_analysis_appendix_C.html#calculatePulsePhase 
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version of the IFGL catalog (jAbdo et alJboiod ) and within 15° of PSR J0034-0534 were modeled 
with power law spectra. The parameters of those point sources > 10° from PSR J0034— 0534 were 
held fixed in the fit. The Galactic diffuse emission was modeled using the glLiem_v02 map cube. 
The extragalactic diffuse and residual instrument background components were modeled jointly 
using the isotropic_iem_v02 template. Both diffuse models are available for download with the 
Fermi STs package. The 7-ray spectrum of PSR J0034— 0534 was modeled as both a power law 
and a simple exponentially cutoff power law, Equation 1 with b = 1, in separate fits: 



dN 
liE 



Nn 



E 



GeYJ 



exp 



(-)' 

\EcJ 



(1) 



Using the likelihood ratio test, the simple exponentially cutoff power law model is preferred 
over a power law at the 4.5(T level. If the emission were to come from very near the stellar surface 
the 7-ray spectrur n could be hyper-exponentia lly cutoff, with b > 1, due to pair attenuation by 
the magnetic field (jPaughertv fc Hardingj|l996l ). Low-altitude emission may be expected given the 
near-alignment of the 7-ray and radio peaks. Assuming a dipolar magnetic field, the predicted 
pair at tenuation cutoff energy is larger than the curvature radiation cutoff energy for nearly all 
MSPs (IHarding. Usov &: Muslimovll2005l ): however, if the field is not dipolar the surface magnetic 
field could be larger and the pair attenuation cutoff energy could be low enough to dominate. The 
spectrum of PSR J0034— 0534 was also fit allowing the b parameter to be free. This fit returned a 
value of b not statistically different from 1 and the 6 = 1 model is still preferred by the likelihood 
ratio test, which is in agreement with the implications of the light curve modeling in §5.1. 

The 7-ray energy spectrum of PSR J0034— 0534, with b = 1, is shown in Figure 2. The plotted 
points in Figure 2 were derived from likelihood fits to each individual energy band in which it was 
assumed the pulsar had a power law spectrum. The energy bands were constructed to be of equal 
size in log space and the last band was chosen to be that which contained the highest energy event, 
6.9 GeV, found consistent with the pulsar position within the 95% containment radius. The best-fit 
parameters are given in Table 2, where the first errors are statistical and the second are systematic. 
The systematic uncertainties were estimated by applying the same fitting procedures described 
above and comparing results using bracketing IRFs where the effective area has been perturbed 
by ±10% at 0.1 GeV, ±5% near 0.5 GeV, and ±20% at 10 GeV with linear extrapolations, in log 
space, between. As a cross check, the 7-ray spectrum of PSR J0034— 0534 was also fit with a binned 
likelihood estimator, ptlike, which computes the photon counts in a point source weighted aperture 
in excess of background counts. The ptlike results agree with the values quoted in Table 2 within 



statistical and systematic errors. Comparison with Table 1 of lAbdo et al.l (l2009al ) shows that the 
spectrum of PSR J0034— 0534 is very typical of the 7-ray MSPs known to date. Table 2 also lists 
the integrated photon flux (F) and energy flux (h) from 0.1 to 100 GeV. 

To search for evidence of modulation at the orbital period events were selected from the data 
described in §3.2 within 5° of the pulsar and with reconstructed energies between 0.1 and 10 GeV. 
The modulation was assumed to be sinusoidal and a maximum likelihood method was used to fit 
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the fraction of the average flux modulated at the orbital period. It might be expected that an 
unpulsed component of emission from, for example, particle acceleration in the wind termination 
shock would show the strongest orbital modulation and thus choosing an off-pulse phase window 
would be best for this search. However, the phase-averaged flux was chosen due to the fact that 
the background estimate shown in Figure 1 does not provide any strong evidence for an unpulsed 
component. The likelihood was constructed by holding the spectral parameters for all sources at 
the phase-averaged fit values and maximizing with respect to the modulated flux and the unknown 
orbital phase of peak emission. There is no evidence for modulation at the orbital period with a 
95% confidence level upper limit on any modulation of 35% of the average flux. 

The Fermi ST gtobssim was used to simulate the region around PSR J0034— 0534 using the fit 
results from the maximum likelihood analysis described in §4.2. The simulation included all point 
sources within 15° of the pulsar and both diffuse backgrounds but did not include the pulsar or a 
model of the 7-ray albedo from the Earth. The simulation start and stop times were matched to 
those in the event file as closely as possible. The same cuts described in §4.1 were applied to the 
simulated data in order to get the background estimates shown in Figure 1. 



DISCUSSION 



5.1. Light Curve Modeling 



In nearly all known 7-ray pulsars with radio counterparts the 7-ray peak(s) lags the radio 
peak(s) by at least 0.05 in phase. The radio profile has tra ditionally been modeled assuming core and 
conal beams centered on the magnetic axis (iRankinll 19931 ) and emitted at low altitude relative to the 
light cylinder radius (Rlc)- Thus, the phase lags between radio and 7-rays have been interpreted as 
indicating different emission altitudes for these tw o wavebands, with the 7-ray emission comi ng from 
outer gaps (OGs) (IRomani Yadigaroglulll995l'l or slot gaps (jMuslimov fc Harding) l2004al ). which 



have a two-pole caustic (TPC) geometry (jPyks &: Rudak 



2OO3I ) ■ reaching much higher altitudes. 



The emission peaks in OG and TPC geometries are due to the formation of caustics at phases where 
relativistic aberration and time-of- flight delays nearly cancel delays from magnetic field curvature 
on trailing field lines (jMorinilll983l ). 



Venter et al. (j2009l ) modeled the light curves of the first eight 7-ray MSPs seen by Fermi 



and observed two distinct subclasses. The radio emission for all eight was fit using a single- 
emission-height conal model. The 7-ray light curves of six MSPs were well modeled using an outer- 
magnetospheric geometry (TPC / OG model), while those of the remaining two were modeled using 
a pair-starved polar cap (PSPC) model in which the 7-ray emission originates from the full open- 
field-line volume in the pulsar magnetosphere, even up to high altitudes. In the TPC / OG case the 
7-ray profile lags the radio, while the radio lags the 7-ray pulse in the PSPC case. In contrast, PSR 
J0034— 0534 is the first MSP for which the radio and 7-ray profiles are observed to be nearly aligned, 
providing strong evidence for co-located emission regions. The term "co-located" is taken to mean 
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that the 7-ray and radio photons are thought to be generated at similar, although not identical, 
locations in the magnetosphere. Fully overlapping emitting regions of the same dimensions will 
lead to identical light curves in these different wavebands, which is not observed. On the contrary, 
the different light curve shapes at different energies seem to imply that the radio emission region is 
smaller, and a subset of the 7-ray emission region in the case of extended, higher-altitude emission. 
The near-alignment of the 7-ray and radio profiles suggests two plausible configurations: either both 
components originate near the polar cap (PC) or both are emitted in the outer magnetosphere. For 
both possibilities the altitude, extent, and degree to which the radio and 7-ray emission regions 
overlap are limited by the shapes of the radio and 7-ray light curves. 

In the first case, relativistic effects tend to smear out the leading peaks while piling up photons 
in the trailing peaks, assuming constant-emissivity annular gaps near the stellar surface extending 
from ~ 0.1 to 0.21HLC- This is contrary to what is observed in the 7-ray light curve. It is also 
difficult to reproduce the symmetric radio peaks assuming emission near the PC, even for emission 
at the stellar surface. If the emission regions are too high in altitude (or their extent is too large) 
the relativistic effects would be boosted (since the corresponding phase shifts due to these effects 
scale as ~ — r/R^c") and the peaks would be even more asymmetric, with a sharp trailing peak. 
Conversely, too little overlap would negate the phase alignment. 

There is however a special solution that is an exception to the rule for the low-altitude models. 
Usually, the peak separation (A) is defined as the difference in phase between the leading and 
trailing peak positions. In certain circumstances, the leading (smeared-out) and trailing (sharp) 
peaks appear "inverted" to an observer at a particular observer angle the angle between the 
rotation axis and the observer's line-of-sight. The trailing peak would be interpreted as the first 
peak, and the leading peak as the second one, since the peak separation between the first and 
second peaks (1 — A in normalized phase) would be smaller than 0.5. If the inclination angle 
between the rotation and magnetic axes (a) is small enough (< 10°), one may find solutions with 
rel ative peak intensit i es an d separations that fit the data quite well. More details are provided 



m 



Venter fc Harding (I2OIOI ). Although this is a valid solution, it is a less likely geometry due to 



the very specific choices of a and C that are needed to obtain a good fit. Unfortunately, the only 
radio polarization measurement of PSR J0034— 0534 to date was inconc lusive and unable to put 
any constraints on the pulsar geometry (jStairs. Thorsett. &: Camildll999l ). 



A co-located outer-magnetospheric origin of both the radio and 7-ray radiation seems more 
likely. If this is the case, PSR J0034— 0534 is the first example of yet another MSP subc lass, 
distinct from the two 7-ray MSP subclasses found bv lVenter et al.l (j2009l ). iManchesterl (j2005l ) has 
suggested that radio emission for young and millisecond pulsars may be generated in the oute r 
magnetosphere, close to the 7-ray emission region, an idea which Weltevrede fc: Johnston (j2008l ) 
explored by comparing the charact eristics of p u lsars with high and low J^sD values in the Parkes 
multi-beam surveys. Additionally, iDyks et al.l (I2OIOI ) have proposed a model in which the radio 
emission for MSPs with symmetric, double-peaked radio features, like what is observed in PSR 
J0034— 0534, arise from curvature radiation from plasma streams with a non-negligible range of 
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emission altitudes. 

Figure 3 shows model 7-ray and radio light curves generated in the context of geometric 'limited 
TPC and OG models', i.e., TPC / OG models with limited extent (along the B-field) of the emission 
regions. The minimum and maximum emission radii for both the radio and 7-ray components are 
constrained by the light curve shapes for each band: excluding lower-altitude emission lowers the 
off-pulse 'shoulder emission' (i.e., emission between consecutive profiles at phases ~ 0.4 — 0.9), and 
'bridge emission' (inter-peak emission at phases ~ 0.0 — 0.3) while giving rise to sharper, more 
separated peaks. On the other hand, including more of the higher-altitude emission boosts the 
leading peak and eventually widens (and to a lesser extent boosts) t he trailing peak. A mor e 



detailed treatment of these 'limited TPC / OG models' can be found in IVenter &: Harding] (|20ld ) 



Given these constraints, the radio and 7-ray light curves are well modeled using a = 30°, 
C = 70°, and a transverse gap width of it; = 0.05 (the fractional angular width starting at the 
PC rim, and normalized to the colatitude of the rim). For the limited TPC models in Figure 3 
the 7-ray emission region extends (in radius) from 0.12Rj;,(7 to O.OR^c' (i-6., starting at the stellar 
surface) while the radio emission region extends from O.GR^c" to O.SHlc- The ranges of emission 
radii for the limited OG models in Figure 3 are the same as those for the TPC model with the 
added caveat that the emission cannot extend below the null charge surface. The 7-ray model 
light curves provide a reasonable fit to the data, although they all somewhat over predict the 
bridge emission as compared to the prediction of the radio curves. The agreement of the observed 
and modeled profiles strongly suggests that both the 7-ray and radio pulses originate in the outer 
magnetosphere with the peaks formed by caustics. It is of interest to note that there is a slight 
excess of events in the 7-ray light curve immediately preceding the first peak which is currently 
insignificant. This feature is reproduced with the TPC model but not with the OG model, since it 
is produced by emission below the null charge surface. With more data, the presence or absence of 
this feature should become statistically better defined and will be useful in further discriminating 
between emission models. 



5.2. Efficiency 

The total luminosity in 7-rays from PSR J0034— 0534 can be calculated using Equation 2: 



ATTfnhd^ 



(2) 



The correction factor depends on the viewing geometry and beaming angle of the pulsar and is 
typica lly ^ 1 for outer- magnetosp heric emission models. Using an approach similar to I Venter et al 
hood ) and IWatters et al.r(|2009l ^ yields fn = 0.74 for TPC models and fn = 0.45 for OG models 



assuming the same geometry used in §5.1. Using the derived energy flux from Table 2 and the 
DM-derived distance from §3.1 yields = 4.7(2.9) ± 3.8(2.3) x lO^^ gj.g g-i foj, ^.j^g TPC(OG) 
model. 
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Once the luminosity is known, the efficiency with which PSR J0034— 0534 converts spin-down 
energy into 7-rays can be calculated using Equation 3: 



(3) 



where Eg/j is the rate at which energy is lost due to the pulsar spinning down. This gives an 
efficiency of r/^ = 0.03(0.02) ± 0.03(0.02) for the TPC(OG) model. Th e derived efficiency for PSR 
J0034— 0534 is on the low end of the efficiencies quoted in Table 1 of lAbdo et al.l (|2009al ). which 
were calculated assu ming fn = 1, and lower than the upper limit of 0.1 estimated for MSPs in 



the GC 47 Tucanae (jAbdo et al.ll2009bl ). The large uncertainties of and rj^ are primarily due to 
the uncertainty of the DM-derived distance in the proper motion. A parallax measurement would 
greatly increase the precision of these values. 



6. CONCLUSIONS 

The MSP PSR J0034-0534 is the ninth 7-ray MSP detected with the Fermi LAT. The 7-ray 
spectral properties of PSR J0034— 0534 are similar to those of other 7-ray MSPs detected by Fermi 
thus far; however, PSR J0034— 0534 is (so far) unique in that it is the only known 7-ray MSP in 
which the radio and 7-ray peaks are very nearly aligned in phase. The near phase alignment can 
be interpreted as providing strong evidence for co-located emission regions. Within the context of 
geometric TPC and OG models, the radio and 7-ray light curves are well modeled by requiring 
that both emission regions be significantly extended in altitude. This implies that both the radio 
and 7-ray emission peaks are a result of caustic formation. A similar radio geometry was required 



to model the Crab 7-ray and radio light curves (iHarding et al.ll2008l ). which are also coincident in 



phase. As the Fermi mission continues and more 7-ray MSPs are detected this new and interesting 
MSP subclass may be further populated. 
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Fig. 1. — The top two panels show the phase- folded light curve of PSR J0034— 0534 for LAT events 
above 1 GeV and above 0.1 GeV within 0.8° of the radio position. 7-ray light curves are shown 
across two rotations with 25 bins per rotation. The dashed horizontal lines correspond to the 
background levels estimated from the simulation described in §4.2. The bottom two panels show 
the Nangay and WSRT radio profiles; the vertical axes are in arbitrary units. 
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Fig. 2. — Phase-averaged 7-ray energy spectrum of PSR J0034— 0534. Plotted points are from 
likelihood fits to individual energy bands where the pulsar is modeled as a power law, solid black 
line is the maximum likelihood model from fitting the full energy range, dashed gray lines are the 
la errors on the model. All sources described in §4.2 were modeled but only the parameters of 
those within 6° were left free in the fits. For each energy band the pulsar was found above the 
background with a test statistic of at least 6, > 2a for two degrees of freedom. 
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Phase (j) 

Fig. 3. — Top: 7-ray data and modeled light curves. Bottom: WSRT 324 MHz radio profile and 
modeled light curves. All modeled light curves were made using a = 30°, C = 70°,andi(; = 0.05. 
The extent of the limited TPC and OG models are given in §5.1. 



Table 1. PSR J0034-0534 Timing Parameters 



Measured Parameters ' 

P (ms) 

P (10-21) 

jjL (mas yr-"*^)^ 



1.8771818845850(2) 

4.966(1) 

31(9) 



Derived Parameters '^ 
d (kpc)'i 

Pcorr(10-21) 
Bsurf (10^ G) 

E5d(1034 erg s-i) 
Blc (10^ G) 



0.53 ± 0.21 
2.63 ± 1.64 
7.12 ± 2.22 
1.57 ± 0.98 
9.89 ± 3.09 



Values in parentheses are rms errors on the last digit. 



TOobbs et al.1 (j2005|). 

■^Derived using Pcorr? whicli accounts for Shklovskii effect. 



^DM-derived using the NE2001 model of lCordes Laziol (|2002l l. 
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Table 2. PSR J0034-0534 7-ray Parameters 
Light Curve Parameters^ 

'^i -0.027 ± 0.008 

FWHMi 0.066 ± 0.019 

51 -0.027 ± 0.008 
02 0.247 ± 0.013 
FWHM2 0.106 ± 0.038 

52 0.011 ± 0.013 
A 0.274 ± 0.015 

Spectral Parameters'^ 

No (10-9 cm-2 s"^ GeV-^) 6.9 ± 1.8 ± 0.2 

r 1.5 ± 0.2 ± 0.1 

Ec (GeV) 1.7 ± 0.6 ± 0.1 

F (10-^ cm-2 s-i) 2.7 ± 0.5 ± 0.4 

/i (10-1^ erg cm-2 s"^) 1.9 ± 0.2 ± 0.1 



^Errors are statistical. Parameters are described in §4.1. 

^First errors are statistical, second are systematic. Parameters are described in §4.2 



